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Abstract 

We argue that the Einstein- Yang-Mills-Higgs theory presents nontrivial solutions with a NUT charge. 
These solutions approach asymptotically the Taub-NUT spacetime and generalize the known dyon black 
hole configurations. The main properties of the solutions and the differences with respect to the asymp- 
totically flat case are discussed. We find that a nonabelian magnetic monopole placed in the field of 
gravitational dyon necessarily acquires an electric field, while the magnetic charge may take arbitrary 
values. 

1 Introduction 

A feature of certain gauge theories is that they admit classical solutions which are interpreted as representing 
magnetic monopoles. For nonabelian gauge fields interacting with a Higgs scalar, there exist even regular 
configurations with a finite mass, as proven by the famous t'Hooft-Polyakov solution 1 . Typically, the 
magnetic monopoles admit also electrically charged generalizations - so called dyons, the Julia-Zee solution 
P] of the SU(2)-Higgs theory possibly being the best known case. These solutions admits also gravitating 
generalizations, both regular and black hole solutions being considered in the literature (see 3 for a general 
review of this topics). In SU(2)-Einstcin- Yang-Mills-Higgs (EYMH) theory, a branch of globally regular 
gravitating dyons emerges smoothly from the corresponding flat space solutions. The nonabelian black 
hole solutions emerge from the globally regular configurations, when a finite regular event horizon radius is 
imposed 01 |S] . These solutions cease to exist beyond some maximal value of the coupling constant a (which 
is proportional to the ratio of the vector meson mass and Planck mass). 

It has been speculated that such configurations might have played an important role in the early stages 
of the evolution of the Universe. Also, various analyses indicate that the monopole and dyon solutions are 
important in quantum theories. 

Since general relativity shares many similarities with gauge theories, one may ask whether Einstein's 
equations present solutions that would be the gravitational analogous of the magnetic monopoles and dyons. 
The first example of such a solution was found in 1963 by Newman, Unti and Tamburino (NUT) [S] [7]. 
This metric has become renowned for being "a counterexample to almost anything" 'S' and represents 
a generalization of the Schwarzschild vacuum solution ^ (see |1(J) for a simple derivation of this metric 
and historical review). It is usually interpreted as describing a gravitational dyon with both ordinary and 
magnetic mass ^. The NUT charge which plays a dual role to ordinary mass, in the same way that electric 
and magnetic charges are dual within Maxwell theory By continuing the NUT solution through its 

horizon one arrives in the Taub universe 0, which may be interpreted as a homogeneous, non-isotropic 
cosmology with the spatial topology S^. 

As discussed by many authors (see e.g. ^lE|)j the presence of magnetic-type mass (the NUT parameter 
n) introduces a " Dirac-string singularity" in the metric (but no curvature singularity) . This can be removed 
by appropriate identifications and changes in the topology of the spacetime manifold, which imply a periodic 
time coordinate. Moreover, the metric is not asymptotically flat in the usual sense although it does obey 
the required fall-off conditions. 

A large number of papers have been written investigating the properties of the gravitational analogs of 
magnetic monopoles ^15.. 16) . the vacuum Taub-NUT solution being generalized in different directions. The 

^ Note that the Taub-NUT spacetime plays also an important role outside general relativity. For example the asymptotic mo- 
tion of monopoles in (super-)Yang-Mills theories corresponds to the geodesic motion in a Euclideanized Taub-NUT background 
1121 . However, these developments are outside the interest of this work. 



corresponding configuration in the Einstein-Maxwell theory has has been found in 1964 by Brill 17 . This 
abelian solution has been generalized for the matter content of the low-energy string theory, a number of 
NUT-charged configurations being exhibited in the literature (see e.g. 18 for a recent example and a large 
set of references). A discussion of the nonabelian counterparts of the Brill solution is presented in |19j . 
These configurations generalize the well known SU(2)-Einstein- Yang-Mills hairy black hole solutions |2U| . 
presenting, as a new feature, a nontrivial electric potential. However, the "no global nonabelian charges" 
results found for asymptotically flat EYM static configurations |2J are still valid in this case, too. 

Here we present arguments for the existence of NUT-charged generalizations of the known EYMH black 
hole solutions Apart from the interesting question of finding the properties of a Yang-Mills- Higgs 

dyon in the field of a gravitational dyon, there are a number of other reasons to consider this type of solu- 
tions. In some supersymmetric theories, closure under duality forces us to consider NUT-charged solutions. 
Furthermore, dual mass solutions play an important role in Euclidean quantum gravity 22 and therefore 
cannot be discarded in spite of their causal pathologies. Also, by considering this type of asymptotics, one 
may hope to attain more general features of gravitating nonabelian dyons. 

The paper is structured as follows: in the next Section we present the general framework and analyse the 
field equations and boundary conditions. In Section 3 we present our numerical results. We conclude with 
Section 4, where our results are summarized. 



2 General framework and equations of motion 
2.1 Action principle 

The action for a gravitating non-Abelian SU{2) gauge field coupled to a triplet Higgs field with vanishing 
Higgs self-coupling is 



S = / ^/^d'x \ - -TUF^^F'-) - -TtiD. Wt-) ) , (1) 



with Newton's constant G. The field strength tensor is given by F^j/ = d^A,, — d^A^ — ie\A^, A^], with 
= d^ — ie[Afj_, ] being the covariant derivative and e the Yang-Mills coupling constant. 
Varying the action Q with respect to g'^" , A^ and <I> we have the field equations 

^ D^i^gFn = ize[$, (2) 



where the stress-energy tensor is 



T^, = 2Tr{F^^F,pg^P - \g^.F^pF'-f} + Tr{]^D^<^D,^ - (3) 
2.2 Metric ansatz and symmetries 

We consider NUT-charged spacetimes whose metric can be written locally in the form 

dr^ 9 
ds^ = -— + P'^(r){de^+sm^edip'^)~N(r)a^(r)(dt + 4nsm^(-^^ (4) 
N[r) 2 

the NUT parameter n being defined as usually in terms of the coefficient appearing in the differential 
dt + 4:n fiiv? {9 / 2) dip. Here 9 and <p> are the standard angles parametrizing an with ranges < 61 < tt, < 
ip < 2-K. 
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Apart from the Killing vector Kq — dt, this line element possesses three more Killing vectors character- 
izing the NUT symmetries 

Ki = sinipde + cos (fi cot 6d^ + 2n cos (fi tan- dt, 

K2 = cos — sin (y9Cot — 2n sin^s tan — 9t, (5) 
Kz = - 2ndt. 

These Killing vectors form a subgroup with the same structure constants that are obeyed by spherically 
symmetric solutions [Ka,Ki,] = eabcKc- 

The nsin^{9/2) term in the metric means that a small loop around the z— axis does not shrink to zero 
at = TT. This singularity can be regarded ar the analogue of a Dirac string in electrodynamics and is 
not related to the usual degeneracies of spherical coordinates on the two-sphere. This problem was first 
encountered in the vacuum NUT metric. One way to deal with this singularity has been proposed by Misner 

His argument holds also independently of the precise functional form of N and a. In this construction, one 
considers one coordinate patch in which the string runs off to infinity along the north axis. A new coordinate 
system can then be found with the string running off to infinity along the south axis with t' = t + Amp, the 
string becoming an artifact resulting from a poor choice of coordinates. It is clear that the t coordinate is 
also periodic with period Sirn and essentially becomes an Euler angle coordinate on S^. Thus an observer 
with (r, 9, if) = const, follows a closed timelike curve. These lines cannot be removed by going to a covering 
space and there are no reasonable spacelike surface. One finds also that surfaces of constant radius have the 
topology of a three-sphere, in which there is a Hopf fibration of the of time over the spatial S'^ (8j. 

Therefore for n different from zero, the metric structure generically shares the same troubles exhibited 
by the vacuum Taub-NUT gravitational field 23 , and the solutions cannot be interpreted properly as black 
holes. 



2.3 Matter fields ansatz 

While the Higgs field is given by the usual form 

$ = (^ra, (6) 

the computation of the appropriate SU{2) connection compatible with the Killing symmetries (01 is a 
more involved task. This can be done by applying the standard rule for calculating the gauge potentials 
for any spacetime group ESj- According to Forgacs and Manton, a gauge field admit a spacetime 
symmetry if the spacetime transformation of the potential can be compensated by a gauge transformation 
CKi^f^ — D^Wi, where C stands for the Lie derivative. 

Taking into account the symmetries of the line element l@J we find the general form 



-^^(^t + Ansir? {^)d(pju{r)T^ + v{r)T'idr + {uj{r)Ti + Cj{r)T2jd9 (7) 

Ja;(r)r2 — J)(r)ri) sin 6* d^p^. 



cos 9t3. 



This gauge connection remains invariant under a residual U{1) gauge symmetry which can be used to set 
— 0. Also, because the variables ui and Q appear completely symmetrically in the EYMH system, the 
two amplitudes must be proportional and we can always set Co = (after a suitable gauge transformation). 
Thus, similar to the n = case, the gauge potential is described by two functions uj{r) and u{r) which we 
shall refer to as magnetic and electric potential, respectively. 



2.4 Field equations and known solutions 

atz, the classical equations of motio: 

dr dt [— - iVP'2 - PP'N') + 2P'{aNPY + — ^) (8) 



Within the above ansatz, the classical equations of motion can be derived from the following reduced action 

,2^3 I 
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'1,,, ,2 Cr(cj2 - 1 + 2nM)2 P2y'2 u^^w^ 1 



where the prime denotes the derivative with respect to the radial coordinate r. 

At this point, we fix the metric gauge by choosing P(r) — \J + , which allows a straightforward 
analysis of the relation with the abelian configurations. 

Dimensionless quantities are obtained by considering the rescalings r — > r/{rje), cf) — s- (prj, n — > n/(rie), 
u rjeu (where rj is the asymptotic magnitude of the Higgs field). As a result, the field equations depend 
only on the coupling constant a = \/ AirGrj. 

The EYMH equations reduce to the following system of five non-linear differential equations 

n^N 1 
rN' = N +-^{-ia'^ -l)-2a^{Nuj'^ + —{uj'^ -l + 2nuf 

rP'' r cr'^JS'^ 

{NaP'^cj)')' = 2auj^(f>, 

( ) = ^(tj^ - 1 + 2nu). 

^ a ' uN P^^ ' 

Two explicit solution of the above equations are well known. The vacuum Taub-NUT one corresponds to 

w(r)=±l, ?i(r) = 0, a(r) = l, (/)(r) = 1, N{r) = 1- '^^^^'^ + f ) _ (iq) 
The U(l) Brill solution 17 has the form 

u;(r) = 0, u{r) =uo+ + ' ^('^) = 1' '^W = 1' = 1 - ^2 + ^2 + 4(^2 ^. „2) " (1^) 

and describes a gravitating abelian dyon with a mass M, electric charge Qe and magnetic charge Qm = 
1 — 2uan, uo being an arbitrary constant, corresponding to the asymptotic value of the electric potential. 

It can be stressed that the Brill solution possesses two, one or zero horizons, according to the values 
of the free parameters Qe, M,u{oo). In the same way as in the case of Reissner-Nordstrom solutions, the 
extremal Brill solution can be defined as the solutions with a degenerate horizon at r = tq. This gives the 
following conditions, fixing M and rg 

2 

To = M , M^+n'-^{Ql + Ql)^0. (12) 

As we will see later, it is convenient to further specify the arbitrary constant u{oo) in such a way the 
u(ro) = 0, this implying 

A'P + n^ ^^■^> 
which fixes Qm and leaves Qe as the only remaining free parameter. In the following we will refer to this 
solution as to the extremal Brill solution. As far as we could see, it is not possible to express M and Qm in 
a closed form depending on (a, n, Qe), but the solution can be constructed numerically. 

2.5 Boundary conditions 

We want the metric (@J to describe a nonsingular, asymptotically NUT spacetime outside an horizon located 
at r = rfi- Here N(rh) = is only a coordinate singularity where all curvature invariants are finite. A 



4 



nonsingular extension across this null surface can be found just as at the event horizon of a black hole. If the 
time is chosen to be periodic, as discussed above, this surface would not be a global event horizon, although 
it would still be an apparent horizon. The regularity assumption implies that all curvature invariants at 
r = r/j are finite. 

The corresponding expansion as r r/i is 



N{r 
a{r 
uj{r 
u{r 
Mr 



Ni{r-rh) + 0{r^rhf, 

<yh + cri(r - r/i) + 0{r - rhf , 

LUh + LUi{r - Th) + 0{r - rhY-, 

ui{r - rh) + U2{r - Vhf + 0{r - rh)^, 

4>h + Mr-rh) + 0{r -rhf, 



(14) 



where = r 



1 

rh 



and 



p2 



U2 



(T/i, ui, LUh, 0h being arbitrary parameters. 
The analysis of the field equations as r 

c, Uq, Qe, (f>i, M 



-^h(ph 



pi 

^h 



CTl 



7lV,,(l-CT,2) 



_!) ^ uiul 



rh 



hVl 



'2ujl(t>h 



2ujIuI 



(15) 



oo gives the following expression in terms of the constants 



N{r) 



2M 


2n^ -a^ 


r 














01 M 


r 





(1 - 2nuf,Y + 



M(2n2 + a202) 



;(r) ~ cc^n/^"^ + . . . , 

Qe n{l - 2nuo) Qe{^r? + o?~^\) 



(16) 



gj,3 



Note that similar to the n ~ Q asymptotically flat case, the magnitude of the electric potential at infinity 
cannot exceed that of the Higgs field, |uo| < 1 ^- The constant M appearing in the asymptotic expansion of 
the metric function N(r) can be interpreted as the total mass of solutions (this can be proven rigurously by 
applying the general formalism proposed in (2H|)- Note that M and n are unrelated on a classical level. 

Also, no purely monopole solution can exist for a nonvanishing NUT charge (i.e. one cannot consistently 
set M = unless w = ±1, in which case the vacuum Taub-NUT solution is recovered). Thus, a nonabelian 
magnetic monopole placed in the field of gravitational dyon necessarily acquires an electric field. 

We close this section by remarking that the definition of the nonabelian charges is less clear for n 7^ 0. 
Although we may still define a 't Hooft field strength tensor, in the absence of a nontrivial two-sphere at 
infinity on which to integrate, the only reasonable definition the nonabelian magnetic and electric charges is 
in terms of the asymptotic behavior of the gauge field. By analogy to the asymptotically flat case, Qe and 
Qra are defined from f'^'' ~ Qe/r"^ and Fg^^ ~ QmSm-O (a similar problem occurs for an U(\) field 1271). 
Thus, since Qm = 1 ^ 2nuQ, the usual quantization relation for the magnetic charge is lost for n ^ 0, which 
is a consequence of the pathological large scale structure of a NUT-charged spacetime. 



■^This depends on the asymptotic structure of the spacetime. 
arbitrary values 1261 . 



For example, in an anti-de Sitter spacetime, uq may take 
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Figure 1. The functions N{r), a{r) and u{r) are plotted for three typical solutions for the same values of (r^, Qe, a). 



3 Numerical results 

Although an analytic or approximate solution appears to be intractable, we present in this section numerical 
arguments that the known EYMH black hole solutions can be extended to include a NUT parameter. 

The equations of motion © have been solved for a large set of the parameters (a, n, Qe, rh), looking 
for solutions interpolating between the asymptotics H14|l and H16|l . NUT-charged solutions are found for any 
n = EYMH dyonic black hole configuration by slowly increasing the parameter n (since the transformation 
n ^ —n leaves the field equations unchanged except for the sign of the electric potential, we consider here 
only positive values of n). As expected, these configurations have many features in common with the n = 
solutions discussed in 0]; they also present new features that we will pointed out in the discussion. Typical 
profiles for the metric functions N{r) and (j{r) and for the electric potential u{r) are presented in Figure 1, 
for a dyonic black hole solution as well as for two NUT-charged solutions. The gauge function w(r) and the 
Higgs scalar ^(r) interpolates monotonically between some constant values on the event horizon and zero 
respectively one at infinity, without presenting any local extremum (see Figure 3). 

The domain of existence of the nonabelian nutty dyons can be determined in the space of parameters. 
If we fix the electric charge Qe of the solution, then there likely exist a volume Vq in the parameter space 
of {a^n,rh) inside which nonabelian solutions exist and on the side of which they become singular and/or 
bifurcate into abelian solution of the type of the Brill solution. For n = the domain of the (a, r^) plane 
where nonabelian solutions exist was determined in |Sj for Qe = and in ^ for Qe 7^ 0. 

The determination of Vq is of course a huge task. In this letter, we will not attempt to determine 
the shape of Vq accurately but rather attempt to sketch it by analyzing the pattern of solutions on some 
generic lines in the space of parameters. For definiteness we set Qe = 0.2 in our numerical analysis, although 
nontrivial solutions have been found also for other values of the electric charge. 

3.1 n varying 

First, we have integrated the system of equations lO with fixed values for a,r/, and Qe and increased the 
NUT charge n. Our values here are a = 1.0, Vh = 0.2 and Qe ~ 0.2 corresponding to a generic values for 
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Figure 2. The values of the parameters M, (j{rh), Nm, Nm, 4>{rh), oj{rh) and it(oo) = uo are shown as a 
function of n for solutions with rn = 0.2, Qe = 0.2 and two different values of the coupling constant a. 

the parameters (the corresponding n = gravitating dyon was constructed in ^). As far as the function 
u{r) is concerned, there exists a main difference between the case n = and n ^ Q. Indeed in the case n = 
this function behave asymptotically like u{r) ^ uq + Qe/r + ©(l/r*^) while in the presence of a NUT charge 
the behaviour is instead u[r) ^ UQ + Qe/r-\-K/r'^, where, as seen from the constant K increases with n. 
Thus, when n becomes large, it becomes more difficult to construct numerical solutions with a good enough 
accuracy ^ , for a given value of the electric charge Qe- 

The effect of increasing n apparently depend strongly of the value a. For a small (typically a < 1) the 
pattern can be summarized by the following points : (i) No local extrema of N(r) are found for small enough 



^To integrate the equations, we used the differential equation solver COLSYS which involves a Newton-Raphson method 
1^ . 
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Figure 3. The metric functions N{r), a{r) and the matter functions w{r), (f){r) and u{r) are shown as a function 
of r for fixed values of [rn, Qe, tl) and three different values of a. The functions N{r) and u{r) of the corresponding 
extremal Brill solution (with uj{r) = 0, h{r) — 1) are also exhibited. 

values of n. For larger n, the function N{r) develops a local maximum and also a local minimum, say 
Nm and N,n at some intermediate, n— depending values of r. For n large enough, we have N^ax > 1- No 
local minimum of N persist for large enough n, the minimum of N(r) (iV^ = 1) being attained as r — > oo. 
(ii) The second metric function a(r) still remains monotonically increasing but the value cr(r/j) diminishes 
when n increases, (iii) The asymptotic value u{oo) also decreases for increasing n. With the values choosen, 
we have u{oo) « 0.189 for n = 0; we find u{oo) = for n « 0.25 and negative values for larger n. 

These effects are illustrated on Figure 2a. On this Figure we have set < n < 2 but we noticed no 
significant change of the behaviour for larger values of n. 

For larger values of a (typically a > 2) the situation is quite different, namely: (i) The function N{r) 
posseses both a local minimum and a local maximum, (ii) The values iy(r/j) and 4>{rh) increase with n and 
approach respectively zero and one, suggesting that the solution approaches an Abelian Brill solution. These 
results are summarize on Figure 2b for a = 2. However, due to numerical difficulties, we could not determine 
properly the value of n where the bifurcation occur. The statement of a bifurcation into a Brill solution is 
confirmed in the next subsection where a is varying. 

Nevertheless, it seems that there are two possible patterns for n oo : for values of a smaller than a 
critical value a, solutions with large values of n seem to occur, while for a > a, the solutions bifurcate into 
a Brill solution (for Qe = 0.2 we find a ~ 1.5). The occurence of these two patterns is reminiscent to the 
case of n = gravitating dyons. 

Note also that, as shown in these plots, the mass parameter M takes negative values for large enough 
values of n. This is not a surprise, since something similar happens already for the U(l) Brill solution Hll|) . 

3.2 a varying 

We now discuss the behaviour of the solutions for a varying a and the other parameters fixed. In absence of 
a NUT charge it is know jSIE] that nonabelian dyonic black hole exist for € [0, y/3 + AQl/2]. For fixed 
Qe and rt and increasing a they bifurcates into an extremal Reissner- Nordstrom solution a,t a ^ Uc- The 
value ac depends of course on rt and Qe- For r/j ^ 1 the value ac ~ 1.4 is found numerically and depends 
weakly on Q^. For th ~ ^?> + AQH2 we have « ^ [i + ^Ql) / {\ + Ql) /2. 

For n > we see (e.g. on Figure 3) that the local maximum characterizing the function N{r) of a nutty 
solution (at least for large enough values of the NUT charge n) progressively disappears in favor of a local 
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Figure 4. The profiles of tlie functions N{r), a{r) and u(r) are represented for Qe ~ 0.2, n = 0.5, a = 1 
and tliree different values of ru ■ 



minimum when a increases. This minimum appears far outside the event liorizon r = rh and becomes 
deeper. In fact, the minimal value N^n approaches zero when a tends to a critical value, say ac{n,Qe,rh)- 
If we denote by r^n, the value of the radial variable where N(rm) = (with > r^) our numerical results 
strongly indicate that the non abelian solution converges into an extremal Brill solution on the interval 
r G [rm,oo) for a — > 

Indeed, the matter functions's profiles u, w, and the metric functions a, N all approaches the profiles 
of the corresponding extremal Brill solution with the same Ofc, Qe, n. This result is illustrated on Figure 3 
for n = 1, rh — 0.2 and Qe ~ 0.2; in this case, we find ac ~ 2.35 but we believe that the result holds for 
generic values of (n, rh, Qe)- The determination of the critical value adn, r/j, Qe) is not aimed in this 
letter. However, it seems that the value Oc depends weakly of n, for example we find Uc ~ 2.22 for n S [3, 4]. 

Nevertheless we can conclude that nutty dyons exist on a finite interval of a and bifurcate into extremal 
Brill solutions for a = ae- 



3.3 Th varying 

In the case of gravitating dyonic black holes with event horizon r/i, the solutions approach the corresponding 
regular gravitating solution on the interval ]r/i,oo[ when the limit rh — > is considered. It is therefore a 
natural question to investigate how nutty-dyons behave in the same limit. 

Considering this problem for a few generic values of (a,n) we reach the conclusion that, in the limit 
rh — > 0, the nutty dyon becomes singular at r = because the value <y{rh) tends to zero. This situation is 
illustrated on Figure 4 where the functions N{r),a{r) and u{r) are plotted for three different values of r^, 
and a = 1, n — 0.5. 

Remarkably, this Figure reveals that the functions cr(r) and u{r) are rather independant of r^ (it is 
also true for w{r),(j){r) which are not represented) while the fimction N{r) indeed involves non trivially 
with rh- Note also that for the metric gauge choice P{r) = Vr^+r?, the area of two-sphere dfl^ = 
P^{r){d9^ + sin^ 9 dip^) does not vanish at r = 0. However, by choosing a Schwarzschild gauge choice 
P{r) = r, a straightforward analysis of the corresponding field equations (which can easily be derived from 
0) implies that it is not possible to take a consistent set of boundary conditions at r = without introducing 
a curvature singularity at that point. Therefore, no globally regular EYMH solutions are found for n 0. 

The determination of the domain of nutty dyons for fixed (a, n, Qe) and increasing the horizon radius r/j 
is very likely an involved problem. Already in the case n — 0, discussed in ^ the numerical analysis reveals 



9 



several (up to three) branches of solutions on some definite intervals of the parameter r^- We believe that 
similar patterns could occur for n > but their analysis is out of the scope of this letter. 

4 Further remarks 

In this work we have analysed the basic properties of gravitating YMH system in the presence of a NUT 
charge. We have found that despite the existence of a number of similarities to the n — case (for example 
the presence of a maximal value of the coupling constant a) , the NUT-charged solutions exhibits some new 
qualitative features. 

The static nature of a n = spherically symmetric gravitating nonabelian solution implies that it can 
only produce a " gravitoelectric" field. There both nonabelian monopole and dyon black hole solutions are 
possible to exist, with a well defined zero event horizon radius limit. For a nozero NUT charge, the existence 
of the cross metric term shows that the solutions have also a " gravitomagnetic" field. The term g^pt 
does not produce an ergoregion but it will induce an effect similar to the dragging of inertial frames j3L)| . 
In this case we have found that only nonabelian dyons are possible to exist and the usual magnetic charge 
quantization relation is lost. The total mass of these solutions may be negative and the configurations do 
not survive in the limit of zero event horizon radius. 

A discussion of possible generalizations of this work should start with the radially excited nutty dyons, 
for which the gauge function w(r) possesses nodes. These configurations are very likely to exist, continuing 
for n > the excited configurations discussed in Also, in our analysis, to simplify the general picture, we 
set the Higgs potential V{(l)) to zero. We expect to find the same qualitative results for a nonvanishing scalar 
potential (at least if the parameters are not to large). It would be a challenge to construct axially symmetric 
NUT-charged dyons (the corresponding n = monopoles are discussed in l^ll). Such dyon solutions would 
present a nonvanishing angular momentum, generalizing the abelian Kerr-Newman-NUT configurations (a 
set of asymptotically flat rotating solutions have been considered recently in |32j ) . 

Similar to the case n — 0, the solutions discussed in this work can also be generalized by including a 
more general matter content. However, we expect that these more general configurations will present the 
same generic properties discussed in this work. This may be important, since there are many indications 
that the NUT charge is an important ingredient in low energy string theory |27 |. conclusion enhanced by 
the discovery of "duality" transformations which relate superficially very different configurations. In many 
situations, if the NUT charge is not included in the study, some symmetries of the system remain unnoticed 
(see e.g. for such an example). Therefore, we may expect the NUT charge to play a crucial role in the 
duality properties of a (supersymmetric-) theory presenting gravitating nonabelian dyons. 

Unfortunately, the pathology of closed timelike curves is not special to the vacuum Taub-NUT solution 
but afflicts all solutions of Einstein equations solutions with "dual" mass in general jSSj. This condition 
emerges only from the asymptotic form of the flelds, and is completely insensitive to the precise details of 
the nature of the source, or the precise nature of the theory of gravity at short distances where general 
relativity may be expected to break down [23) . This acausal behavior precludes the nutty dyons solutions 
discussed in this paper having a role classically and implies a number of pathological properties of these 
conflgurations. 

Nevertheless, there are various features suggesting that the Euclidean version of NUT-charged solutions 
play an important role in quantum gravity |22| . For example, the entropy of such solutions generically do 
not obey the simple " quarter-area law" . As usual, a positive-definite metric is found by considering in 
the analytical continuation t it, n —> i n, which gives P^(r) = — r? . In this case, the absence of conical 
singularities at the root of the function N(r) imposes a periodicity in the Euclidean time coordinate 

47r 

, , y (17) 

7V'(r/i)o-(r/i) 

which should be equal with the one to remove the Dirac string fi — Sim. In the usual approach, the solution's 
parameters must be restricted such that the fixed point set of the Killing vector dt is regular at the radial 
position r — rfi- We find in this way two types of regular solutions, "bolts" (with arbitrary — ri, > n) or 
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" nuts" (r/i = n) , depending on whether the fixed point set is of dimension two or zero (see |13| for a discussion 
of these sohitions in the vacuum case and |34| for a recent generahzation with anti-de Sitter asymptotics). 

We expect that the EucHdean nutty dyons wiU present some new features as compared to the Lorentzian 
counterparts. For example, globally regular solutions may exist in this case, since r = corresponds to the 
origin of the coordinate system (note also that the SU(2) Yang-Mills system is known to present self-dual 
solutions in the background of a vacuum Taub-NUT instanton j33|)- In the absence of closed form solu- 
tions, the properties of these non-self dual EYMH solutions cannot be predicted directly from those of the 
Lorentzian configurations. However, similar to the Lorentzian case, they can be studied in a systematic way, 
by using both analytical and numerical arguments. For example, the magnitude of the electric potential at 
infinity of the Euclidean solutions, is not restricted. Also, the condition /3 — Snn implies N'{rh)<j{rfi) — 2n 
and introduces a supplementary constraint on the matter functions as r — s- r/j. A study of such solutions 
may be important in a quantum gravity context. 
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